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Summary. This talkQ will review selected topics in rapidity gap physics. In par- 
ticular I will discuss diffractive jet production and the possibility of searching for 
the higgs boson using diffraction at the LHC; the dipole picture of diffraction and 
saturation; and those processes where a large momentum is transferred across the 
rapidity gap, for which there has been recent progress both experimentally and 
theoretically. 

1 Introduction 

Over the past 10 years, due in no small part to the quality and extent of 
data collected at the HERA and Tevatron colliders, the field of rapidity gap 
physics has flourished. As a result, in a review talk like this I cannot hope to 
cover anything other than a few topics, chosen to reflect my personal bias. 

The next section will focus on the hard diffractive production of jets and 
higgs bosons, and will draw on data collected at both HERA and the Teva- 
tron. In Section 3, I discuss the dipole model of diffraction and the evidence 
for saturation. In Section 4, I turn to rapidity gaps with a large momentum 
transfer across the gap. These are rarer but rather clean processes, and recent 
high quality data on vector meson production has allowed comparison with 
theory, which I discuss. 

2 Hard diffraction 



I will follow the conventional language and classify hard diffractive processes 
as shown in Figure |l|. Q is some hard scale (e.g. jet transverse momentum, 
W mass etc.) characteristic of the system X. In "double pomeron exchange" , 
which has been measured at the Tevatron jy, the protons remain intact, 
losing only a small fraction of their initial energy, so that the system X is 
produced centrally. In single diffraction, only one proton remains intact and 
fast, and the system X is distant from it in rapidity. 

1 Talk presented at the 14th Topical Conference on Hadron Collider Physics, HCP 
2002, Karlsruhe, 29 September-4 October 2002. 
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"[iQiible pomeion exchange" "single diffraction" 



Fig. 1. Hard diffraction 




Fig. 2. Regge factorization. (Figure from W) 



At HERA, single diffraction of a virtual photon (virtuality Q 2 ) has allowed 
experimenters to probe the partonic structure of the diffractive exchange [S|| . 
Appealing to regge theory, one can attempt to write the cross-section as a 
product of a pomeron flux factor fp/ p and a pomeron parton density function, 
as shown in Figure |^. The experimenters parameterise the flux factor as 

fp/ P (xp,t) = aa(t) _ 1 (1) 
x p 

where xp is the fraction of the incoming proton's energy carried by the 
pomeron, t is the momentum transfer to the scattered proton, B is the diffrac- 
tive slope and a(t) is the pomeron trajectory. The experimenters are able to 
fit all their data on the diffractive structure function using a parameterisa- 
tion of this form after evolving the parton density functions using the NLO 
evolution equations. HI finds a pomeron intercept ap = 1.17 ± 0.02 and no 
need for secondary regge exchanges for xp < 0.01 ||. HI has now extracted 
the pomeron quark and gluon density functions with an estimate of the error, 
see Fig. [| Note that it is the gluon at large z = x/xp which is least well 
constrained. 

If the notion of a universal pomeron parton density function is to be 
tested, one needs to take the partons as measured in diffractive DIS, and use 
them to predict the rates for other processes. At HERA this programme has 
been carried out quite extensively and on the whole there is good agreement 
[|J. However, there are indications that the DIS partons do tend to lead 
to an overestimate of the data collected off real photons. Good agreement 
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Fig. 3. Pomeron partem density functions extracted by HI [^) 



can be arranged if one is prepared to accept an overall renormalisation by a 
factor of about 0.6 ||. Strictly speaking, the need for such a renormalisation 
violates universality. This should not come as a surprise; we already knew 
that universality should not hold across the board in diffraction M. What 
would be of interest is if the violation can be understood. Simple (eikonal) 
models predict that rapidity gaps will be filled in by secondary interactions in 
those processes where the incoming beam particles have structure, which is 
the case in hadron-hadron interactions and photon-hadron interactions with 
an on-shell photon. These simple models also predict that the filling in of 
gaps can be approximated by an overall multiplicative factor which is weakly 
process dependent; depending primarily on the overall centre-of-mass energy 

In this way we can understand a gap survival factor of ~ 0.6 at HERA. 
The eikonal models also predict a survival factor of around 0.1 at the Tevatron 
(the high centre-of-mass energy being the main reason for the reduction since 
it liberates more low- a; partons). 

The burning question is therefore: "How does a gap survival factor of 
0.1, in conjunction with the latest HI parton density functions, stand up to 
the Tevatron data?". At first sight, the answer is "very badly". In Table |], 
we show the original calculations of which are parton level and do not 
contain any gap survival factor. Comparison is to CDF and DO data available 
at the time (some of which were preliminary) . The references shown in the 
table are to the final published papers. The pomeron parton densities do not 
now agree with the most recent HERA data, but they are not so far out to 
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account for the obvious problems. Even with a gap survival of ~ 0.1, things 
look bleak. 



Table 1. Problems at the Tevatron? 





Experiment 


Theory 


Theory /Exp. 


CDF W @ 


1.15 ±0.51 ±0.20% 


7% 


6±3 


CDF RG dijet @ 


0.75 ± 0.05 ± 0.09% 


16% 


22 ± 3 


CDF pot dijet @ 


0.109 ± 0.003 ±0.016% 


4% 


34 ±5 


DO RG dijet fl| 


0.67 ± 0.05% 


12% 


18 ±1 


CDF heavy quark jjj 


0.18 ±0.03% 


30% 


167 ± 28 


CDF double pomeron exch. fi] 


13.6 ±2.8 ±2 nb 


3713 nb 


273 ± 69 
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Fig. 4. Comparison of the diffractive dijet rate to CDF data. (Figure from H 



However, the problem may not be so bad. With the new HI partons 
the rate for diffractive dijet production at the Tevatron agrees well with the 
HERA prediction with a gap survival factor of 0.1. The Tevatron data on 
diffractive dijet production now go beyond the total rate, as can be seen in 
Fig. U {Ffj is a ratio of diffractive to non-diffractive cross-sections and (3 is the 
momentum fraction of the parton coming from the pomeron which is involved 
in the hard subprocess). The solid red curve shows the prediction based on 
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the latest HI partons - it agrees in shape with the CDF data, especially once 
one realises that the large j3 region is the large z region in Fig. || and that 
the dijet process is gluon dominated. So much for dijet production in single 
diffraction. What about the double pomeron exchange process (producing 
central dijets), which is out by two orders of magnitude according to Table 

Remarkably even here things seem not to be too bad. Analysis shows 
that there are large hadronisation corrections to the parton level results of 
Table [l] which account for a suppression of the theory by a factor of 4 |H . 
This suppression comes about because CDF used a small cone R = 0.7 to 
define their central jets, simultaneously with a low Et cut of 7 GeV. At such 
low scales the jets are broad and one loses about 2 GeV per jet. In a 
HERWIG monte carlo simulation of the parton showering and hadronisation 
[ fl6| was included thereby allowing us to estimate the size of the corrections. 
In addition, the overall rate is sensitive to the pomeron intercept. A higher 
intercept leads to a lower cross-section at the Tevatron if the normalisation 
is fixed at HERA since HERA probes lower values of xp. The authors of M 
used a soft pomeron intercept which we now know to be inappropriate; this 
gains another factor of 2. The remaining difference is down to the parton 
densities in the pomeron. Table | summarises the results of (l5) where the 
bottom line shows a difference of a factor 10 between theory and experiment 
which is in accord with gap survival estimates. Note that there is good reason 
to expect significant "contamination" from secondary exchanges (IR) in the 
region of the Tevatron data. 



Table 2. Rates for double pomeron exchange: comparison of theory and experi- 
ment. Theory calculations performed using HI fit 3 (LO) partons and the default 
secondary exchange of pM 



Regge exchange 


Parton Level (nb) 


Hadron Level (nb) 


P 


1175 


339 


R 


241 


58 


P + R 


1416 


397 


CDF data 




43.6 ±4.4 ± 21.6 



2.1 Central higgs production 

It has been suggested that the double pomeron exchange process could be 
utilized at the LHC to produce new particles, in particular the higgs boson 
[ pT[p^ | . If the higgs is produced exclusively, as shown in Fig. ||, then one could 
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Fig. 5. Exclusive higgs production. (Figure from ]l8|]) 

reconstruct its mass quite accurately (to within 1 GeV (l9)) by tagging the 
outgoing protons. Moreover, the exclusive nature of the central system leads 
to a significant suppression of QCD backgrounds, so that one could utilise the 
bb decay of the higgs with a much better S/ B compared to the non-diffractive 
production mechanism. Khoze, Martin and Ryskin (KMS) calculate the rate 
for the diagram in Fig. [| and find 3 fb for a 115 GeV higgs decaying to bb at 
the LHC, which should be sufficiently large to permit a good measurement. 
It is possible to test the reliability of this estimate by performing the cor- 
responding calculation for exclusive central dijet production (i.e. replace the 
higgs by a dijet pair). KMS predicted a rate of around 1 nb at the Tevatron, 
which is to be compared to the CDF upper limit of 4 nb |Q, i.e. the process 
has not been seen in Run I data. However, with the increase in data from 
Run II, and the fact that DO now has roman pot detectors in both forward 
and backward directions, it should be possible to check the KMS calculation. 

3 Dipole models 



Dipole 



Fig. 6. Diffraction of a colour dipole 



Rapidity Gap Physics 7 



Complementary to the regge picture of diffraction in photon induced re- 
actions is the dipole picture, to which we now turn our attention. In the 
proton rest frame, the incoming photon converts into a qq pair a long dis- 
tance upstream of the proton. In the diffractive limit, the quarks are highly 
energetic and travel along straight lines through the proton, picking up a 
non-Abelian phase factor, before eventually forming the diffracted system X 
way downstream of the proton, see Fig. 0. Consequently, one can write down 
expressions for a variety of diffractive process, e.g. for the total jp cross- 
section at high-energy we only need the imaginary part of the forward elastic 
scattering amplitude, i.e. 

a? p L = J dzd 2 r\^ L (z,r)\ 2 a( S ,r,z) (2) 

where a is the cross-section for scattering a colour dipole of transverse size r 
and energy fraction z off a proton. It is universal in that the same dipole cross- 
section should appear in other processes, such as diffractive vector meson 
production, where one simply replaces the outgoing photon wavefunction with 
the meson wavefunction |2(|. There is clearly a lot of physics in the dipole 
cross-section, including the QCD evolution of the original dipole. Apart from 
the total jp cross-section (and hence the structure functions, F%, Fl and 
F£) and vector meson production, the dipole formalism has been used to 

compare to data on inclusive diffractive DIS (kP^), deeply virtual Compton 
scattering and shadowing off nuclei. 

In the dipole model of Golec-Biernat & Wiisthoff pj[ , the dipole cross- 
section was parameterized as 



where 



and 



R °( x ) = 77-T72 \ 

GeV V^o 



Q 2 + 4m 2 . 



W 2 ' 

Ro(x) is called the saturation radius, since for larger r the cross-section flat- 
tens off. Since the saturation radius moves to smaller r as x decreases this 
model naturally tames the powerlike behaviour of the total cross-section as 
one moves to smaller x. For small enough r, the cross-section goes like r 2 
which generates Bjorken scaling. The striking agreement of the model with 
data i*2 and on F®^ originally led to the idea that HERA was already prob- 
ing the non- linear dynamics of saturation. Subsequent, more detailed studies 
using the latest data, have revealed that it is not possible to fit the data with 
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a pure power, i.e. A, which gets tamed by saturation effects. It is necessary 
to replace the exponent with the gluon density p2| |, i.e. 

The gluon density itself becomes less steep as Q 2 falls and so the exponenta- 
tion is less important. This means that the regime of large corrections arising 
from the non-linear dynamics is pushed beyond the HERA region, i.e. to 
smaller x. The need to go beyond the original model of Golec-Biernat & 
Wiisthoff is illustrated in Fig. fj] where the effective slope (~ x _A ) is shown as 
a function of Q 2 . In the original model, this slope asymptotes to just below 
0.3 and deviations at lower Q 2 are wholly attributed to saturation dynamics. 
In the new model, there is no flattening at high Q 2 . 




Fig. 7. The effective slope of the low x structure function. Solid curve is the new 
dipole model, whilst the dashed curve is that from the original Golec-Biernat & 
Wiisthoff model. (Figure from Q) 



A quite different approach can be found in ]23| . Here the dipole cross- 
section is written as a sum of two terms each of which can be thought of as 
arising from a pure regge pole (of intercepts 1.06 and 1.4 respectively). This 
"two pomeron" model has no saturation dynamics at all. The reduction of 
the effective slope at low Q 2 arises because of the dominance of the pomeron 
with lower intercept in that region (the reverse occuring at high Q 2 ). This 
model is also able to describe the available data, including the diffractive 



structure function p4| and deeply virtual Compton scattering p5 
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Before leaving dipoles, I should say a few things about the latest theo- 
retical progress in the physics of saturation. The use of non-linear perturba- 
tive QCD dynamics to control the growth of low- a; cross-sections has a long 
history, dating back to the "GLR equation" of the early 1980's pfl. More 
recently, the Balitsky-Kovchegov equation has been developed to describe 
the non-linear evolution of the S-matrix for scattering a colour dipole off a 
hadronic target [^7j . Underpinning all of this is the colour-glass-dynamics of 
[ ^8| . Formulated as an effective field theory (analogous to that of glasses in 
condensed matter physics) , the colour glass dynamics describes the quantum 
evolution of soft gluons in a classical background colour field. It reduces to 
the BFKL equation in the approximation of a dilute background, and to the 
Balitsky-Kovchegov equation in the large N c limit. 

4 Rapidity gaps at high-i 

So far, all the processes we have looked at are close to t = 0, i.e. the outgoing 
proton(s) do not receive a large transverse momentum. Let us now focus on 
the case the —t ^> A 2 . In this case, the incoming proton will typically be 
shattered. It is thought that the largeness of the momentum transfer will 
allow us to utilise QCD perturbation theory and hence to test the relevance 
of BFKL dynamics in these processes. 

4.1 Vector mesons 




Fig. 8. High pr vector meson production. (Figure from [p0[) 

In Fig. |^ we show the diffractive production of a high px vector meson, V. 
The large momentum transfer across the gap, —t re pj~, almost always breaks 
up the proton and leads to a jet in the resultant debris. The non-perturbative 
dynamics is factorised into either the proton parton density functions or into 
the meson lightcone wavefunction, and QCD can be used to compute the 
dynamics of the exchange. The fact that xW 2 3> —t ensures that there is 
still a large rapidity gap between the proton dissociation products and the 
vector meson. 
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In leading order, BFKL predicts that the hard subprocess cross-section 
for scattering off quarks (gluons differ only by a colour factor) should go like 

m 



da(jq^Vq) a 4 s e Szla2 ( In 2 



T 



dt t 4 z3/ 2 * \ 112zC(3) 

where 

-t , 3a s , xW 2 
and z = In ■ 



ex P — TToTTToT . ( 5 ) 



Q 2 + my 27r s 



10 3 




Fig. 9. Comparison of data |32j with theory for high pr p production. The solid 
line is LO BFKL and the dotted line two-gluon exchange. (Figure from fjif) 

Comparison of the leading order BFKL calculation with the ZEUS data is 
shown in Fig. || |3jJ. The calculation is a fit to the data with so = f3my — jt 
and a s being the free parameters. The fit shown corresponds to j3 = 0, 7 = 1 
and a s = 0.20. Good agreement is also found (with the same parameters) 
with the data on the (j> and J/ ip mesons |33| . Note that it is not possible to 
get good agreement in the approximation that only two-gluons are exchanged 
between the diquark system and the struck parton. 

The above curves were computed assuming a very simple form for the me- 
son wavefunction. In particular, it is assumed that the quarks share equally 
the momentum of the meson. Relativistic corrections to this simple approxi- 
mation have been considered [^9| and do not appear to spoil the good agree- 
ment j30|. Inclusion of relativistic corrections also allows one to quantify the 
degree to which the helicity of the meson differs from that of the photon, and 
comparison with data is once again encouraging ]30| ] . 

In the future, data will become available on high-t photon production. 
This is something to look forward to, since it avoids the uncertainty associated 
with the production of the vector meson. 
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4.2 Gaps between jets 



A7J 




Fig. 10. Dijet production with a rapidity gap 



As well as vector meson production, one can look for rapidity gaps between 
jets in photon-hadron and hadron-hadron collisions |$4j. The typical final 
state topology is shown in Fig. where two jets are produced far apart in 
rapidity and there is a gap between the jets. Early measurements at HERA 
f[35| and the Tevatron |3(|[37j have been compared to theory and leading order 



BFKL does fine J38 39 . However, conclusive statements are hard to make 



either because the gap is not large enough (i.e. the excess over non-BFKL 
QCD is small) or hadronisation corrections are large pq. 
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Fig. 11. Recent HI results on gaps between jets shown as a function of rapidity 
gap and transverse energy between the jets. (Figure from |40| ) 



There has been significant recent progress, both experimentally and the- 
oretically in this area which has to some extent shifted interest away from 
BFKL dynamics and into the domain of jet energy flows. HI has focussed on 
the definition of the gap. They use the kr cluster algorithm to put all hadrons 
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into jets, after which they select the two highest pr jets. The summed trans- 
verse energy between these two jets (Z?f, ap ) is then used to define a gap. As 
can be seen in Fig. O, for low i?f, ap , one can really speak of a rapidity gap 
and a very clear excess is seen in the data over the standard Monte Carlos, 
whilst at larger values of £?f, ap the enhancement is less pronounced . By 
defining their gaps this way, HI has reduced its sensitivity to soft gluon ra- 
diation by effectively cleaning up the edge of the gap in a way which makes 
possible direct comparison with future theoretical calculations. 

On the theoretical side, Dasgupta & Salam have recently pointed out that 
there is a previously unconsidered mechanism which ought to be considered 
when considered interjet energy flows p~ij| . In particular, they have discovered 
a class of "non-global" logarithms which ought to be summed at the single 
logarithm level. 

5 Summary 

The key conclusions of this talk can be summarised as follows: 

• Regge factorisation and QCD evolution work well to describe diffractive 
deep inelastic scattering at HERA 

• HERA partons may well be useful at the Tevatron provided one accounts 
for gap survival at a level of around 10%. Tevatron measurements at lower 
xp would help. 

• There is a need to measure exclusive central dijet production at Tevatron 
Run II in order to constrain diffractive higgs cross-sections at the LHC. 

• Position on saturation is not clear. Ideally one would like to go to lower 
x. 

• High-i vector meson production is well described by leading order BFKL. 
Acknowledgements 

Special thanks to Brian Cox, Rikard Enberg and Paul Newman for their 
help in preparing this talk. Thanks also to the conference organisers for their 
invitation to participate. 

References 

1. T. Affolder et al, CDF Collaboration: Phys. Rev. Lett. 85, 4215 (2000). 

2. C. Adloffet al, HI Collaboration: Z. Phys. C 76, 613 (1997). 

3. J. Breitweg et al, ZEUS Collaboration, Eur. Phys. J. C 6, 43 (1999). 

4. P. Thompson: 'Jets and Et flows in diffraction at HERA', talk presented at DIS 
2002, April-May 2002, Acta Phys. Polon B 33, 3213 (2002). 

5. F.-P. Schilling: 'Measurement and NLO DGLAP QCD Interpretation of Diffrac- 
tive DIS at HERA', talk presented at DIS 2002, April-May 2002, Acta Phys. 
Polon B 33, 3511 (2002). 



Rapidity Gap Physics 13 



6. S. Schatzel, 'Diffractive final states at HERA', talk presented at the Workshop 
on Low x Physics, Antwerpen, September 2002. 

7. J.C. Collins: 'Factorization in hard diffraction', talk presented at the Ringberg 
Workshop, J. Phys. G 28, 1069 (2002), and references therein. 

8. E. Gotsman, E. Levin, U. Maor: Phys. Rev. D 60, 094011 (1999); A. Kaidalov, 
V.A. Khoze, A.D. Martin, M.G. Ryskin: Eur. Phys. J. C 21, 521 (2001). 

9. L. Alvero, J.C. Collins, J. Teron, J.J. Whitmore: Phys. Rev. D 59, 074022 
(1999); talk presented at DIS 1998, April 1998; L. Alvero, J.C. Collins, J.J. 
Whitmore: Nucl. Phys. Proc. Suppl. 79, 382 (1999). 

10. F. Abe et al, CDF Collaboration: Phys. Rev. Lett. 78, 2698 (1997). 

11. F. Abe et al, CDF Collaboration: Phys. Rev. Lett. 79, 2636 (1997). 

12. T. Affolder et al, CDF Collaboration: Phys. Rev. Lett. 84, 5043 (2000); T. 
Acosta et al, CDF Collaboration: Phys. Rev. Lett. 88, 151802 (2002). 

13. B. Abbott et al, DO Collaboration: Phys. Lett. B 531, 52 (2002). 

14. T. Affolder et al, CDF Collaboration: Phys. Rev. Lett. 84, 232 (2000). 

15. R. Appleby, J.R. Forshaw: Phys. Lett. B 541, 108 (2002). 

16. B.E. Cox, J.R. Forshaw: Comput. Phys. Commun. 144, 104 (2002); 
www.pomwig.com. 

17. A. Schafer, O. Nachtmann, R. Schopf: Phys. Lett. B 249, 331 (1990); A. Bialas, 
PV. Landshoff, Phys. Lett. B 256, 540 (1991); Yu. L. Dokshitzer, V.A. Khoze, 
T. Sjostrand: Phys. Lett. B 274, 116 (1992); J.D. Bjorken: Phys. Rev. D 47, 
101 (1993). 

18. V.A. Khoze, A.D. Martin, M. G. Ryskin: 'Diffractive Higgs Production: Myths 



and Reality', |hep-ph/0207313| Eur. Phys. J. C 14, 525 (2000). 

19. A. de Roeck, V.A. Khoze, A.D. Martin, R. Orava, M.G. Ryskin: Eur. Phys. J. 
C 25, 391 (2002). 

20. N.N. Nikolaev, B.G. Zakharov: Z. Phys. C 53, 331 (1992); B.Z. Kopeliovich, J. 
Nemchik, N.N. Nikolaev, B.G. Zakharov: Phys. Lett. B 309, 179 (1993); N.N. 
Nikolaev, B.G. Zakharov: Z. Phys. C 64, 631 (1994); J. Nemchik, N.N. Nikolaev, 
E. Predazzi, B.G. Zakharov: Z. Phys. C 75, 71 (1997). 

21. K. Golec-Biernat, M. Wiisthoff: Phys. Rev. D 60, 114023 (1999); Eur. Phys. J. 
C 60, 313 (2001). 

22. J. Bartels, K. Golec-Biernat, H. Kowalski: Phys. Rev. D 66, 014001 (2002). 

23. J.R. Forshaw, G. Kerley, G. Shaw: Phys. Rev. D 60, 074012 (1999). 

24. J.R. Forshaw, G. Kerley, G. Shaw: Nucl. Phys. A. 675, 80c (2000). 

25. M.F. McDermott, R. Sandapen, G. Shaw: Eur. Phys. J. C 22, 655 (2002). 

26. L.V. Gribov, E.M. Levin, M.G. Ryskin: Phys. Rep. 100, 1 (1983). 

27. Ia. Balitsky: Nucl. Phys. B 463, 99 (1996); Yu. Kovchegov: Phys. Rev. D 60, 
034008 (2000); ibid D 61, 074018 (2000). 

28. see E. Iancu: ' The Colour Glass Condensate', plenary talk at Quark Matter 
2002, July 2002, |hep-ph/0210236| , and references therein. 

29. D.Yu. Ivanov, R. Kirschner, A. Schafer, L. Szymanowski: Phys. Lett. B 478, 
101 (2000); erratum ibid B 498, 295 (2001). 

30. R. Enberg, L. Motyka, G. Poludniowski: talk presented at DIS 2002, Cracow, 
Poland, April-May 2002, Acta Phys. Polon B 33, 3511 (2002). 

31. J.R. Forshaw, M. Ryskin: Z. Phys. C 68, 137 (1995). 

32. S. Chekanov et al, ZEUS Collaboration: DESY-02-072, to be published in Eur. 
Phys. J. C. 

33. J.R. Forshaw, G. Poludniowski: hep-ph/0107065 , to be published in Eur. Phys. 
J. C. 



14 Jeffrey R. Forshaw 



34. A.H. Mueller, W.-K. Tang: Phys. Lett. B 284, 123 (1992). 

35. M. Derrick et al, ZEUS collaboration: Phys. Lett. B 369, 55 (1996). 

36. F. Abe et al, CDF Collaboration: Phys. Rev. Lett. 74, 855 (1995); ibid 80, 
1156 (1998); ibid 81, 5278 (1998). 

37. S. Abachi et al, DO Collaboration: Phys. Rev. Lett. 76, 734 (1996); ibid 84, 
5722 (2000). 

38. B.E. Cox, J.R. Forshaw, L. Lonnblad: JHEP 9910:023 (1999). 

39. R. Enberg, G. Ingelman, L. Motyka: Phys. Lett. B 524, 273 (2002). 

40. C. Adloffet al, HI collaboration: Eur. Phys. J. C 24, 517 (2002). 

41. M. Dasgupta, G. Salam: JHEP 0203:017 (2002). 



